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In this article the extended Bose-Hubbard model describing ultra-cold atoms confined in a shallow, 
one-dimensional optical lattice is introduced and studied by the exact diagonalization approach. All 
parameters of the model are related to the only relevant parameter controlled experimentally - the 
depth of the optical potential. Changes in a shape of the insulating lobe in the phase diagram of 
the system are explored and the value of the critical tunneling for which the system undergoes the 
phase transition (from the insulating to the snperfluid phase) is predicted. It is shown that the 
value of critical tunneling is substantially affected by the presence of the tunnelings to distant sites 
of the optical lattice. The results may have some importance in upcoming experiments on quantum 
quench through phase transition points. 


I. INTRODUCTION 

The celebrated Hubbard model was introduced over 
fifty years ago to give an intuitive explanation of the tran¬ 
sition from the superconducting to the insulating phase 
for not fully filled conducting bands [1]. The model re¬ 
vealed that on the many-body level the quantum phase 
transitions are induced not only by a mobility of single 
particles but also by mutual interactions between them. 
Due to amazing experimental progress in controlling of 
ultra-cold gases in optical lattices, this old fashioned the¬ 
oretical toy-model is undergoing a renaissance since it 
provides a realistic description of the real quantum sys¬ 
tems m]. Many different extensions of the standard 
Hubbard model came into a play due to a subtle analysis 
of interactions between particles. When interactions are 
strong enough they can induce inter-site and inter-band 
couplings and may lead to many interesting phenomena 
[5HI3]. This path was intensively explored by many au¬ 
thors and recently it was summarized in a broad review 
|I4j . Nevertheless, there are still open questions on the 
properties of ultra-cold gases confined in a shallow op¬ 
tical lattice. One suspects that in such a case the va¬ 
lidity of the Hubbard-like Hamiltonian can be question¬ 
able and some other methods should be adopted. One 
possibility is to exploit some methods working directly 
in the configuration space of confined particles. This 
idea was adopted recently with so-called hybrid quan¬ 
tum Monte Carlo method and used to study the phase 
transition from the Mott-Insulator (MI) to the superfluid 
(SF) phase [15] . 

In this article the old-fashioned Hubbard-like approach 
is applied to this problem. First, we start from the gen¬ 
eral Hamiltonian of the system and we rewrite it in the 
form of the exact multi-orbital Hubbard like model. Then 
we show that even for very shallow lattices, for the lowest 
commensurate filling of the lattice, p = 1, in the vicinity 
of the phase transition the single-particle tunneling to 
the next-nearest-neighbor (NNN) sites has to be taken 
into account. Moreover, we show that NNN tunneling 
is a dominant correction to the standard Bose-Hubbard 


(BH) model and in the first approximation it is sufficient 
to neglect other corrections. In this way we obtain a quite 
simple extended Bose-Hubbard (BH) model and we study 
its properties in the vicinity of the phase transition with 
the exact diagonalization method. 

The paper is organized as follows. In Sec. jH] the ex¬ 
tended BH model with tunnelings to the NNN is intro¬ 
duced. In Sec. |HI| the method of the exact diagonaliza¬ 
tion of the Hamiltonian is explained and an extrapolation 
of the results to the thermodynamic limit is clarified. In 
this way the phase diagram of the system for different 
lattice depths is calculated and the critical tunneling is 
localized. Finally, the conclusions are presented in Sec. 
IlYl 


II. THE MODEL 

The starting point of the derivation of the model is a 
general Hamiltonian describing bosons of mass m con¬ 
fined in an optical lattice potential and interacting mu¬ 
tually via short range delta-like potential. We assume 
that the dynamics of the system is completely frozen in 
directions perpendicular to the direction of the optical 
lattice. Experimentally it can be obtained by applying 
very strong harmonic confinements in these directions. 
Then, excitations to higher states of harmonic confine¬ 
ment are strongly suppressed and, in consequence, parti¬ 
cles occupy only the ground-state of the harmonic oscil¬ 
lator. Effectively, the second-quantized Hamiltonian of 
the system in the direction of the optical lattice has a 
form: 

H = J dx^\x)HQ'^{x)^ j dx'F^(a;)'I'^(x)d'(a;)'I'(a;), 

( 1 ) 

where + V sm^{kx) is a single-particle 

part of the Hamiltonian. Lattice intensity V and the 
wavevector k are controlled by the lasers forming an op¬ 
tical lattice. The effective one-dimensional contact in¬ 
teractions are controlled by the coupling constant g US]. 
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The bosonic field operator annihilates a particle at 
point X and fulfills the standard commutation relations 
for bosonic fields = 5{x — x’). For conve¬ 

nience the whole analysis will be performed in the natural 
units of the problem, i.e. we measure energies in recoil 
energies /2m, lengths in 1/k, etc. In this way 

the coupling constant g is dimensionless and it measures 
an effective strength of the interaction between particles. 

The standard route to obtain the effective Hubbard¬ 
like model describing the system is to expand the field 
operator in the basis of the maximally localized Wan- 
nier functions 'W[°‘\x). The Wannier functions are num¬ 
bered with two discrete indices a and i corresponding 
to the Bloch band and lattice site of the periodic poten¬ 
tial, respectively. In the case studied, the functions can 
be found straightforwardly by solving the single-particle 
Schrodingier equation, which has a form of well-known 
Mathieu equation. The expansion of the field operator 
in the basis of Wannier functions has a form 
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^(x) =^W^°\x)aai, ( 2 ) 

i,Oi 

where a^i annihilates a boson in the single particle state 
described by the wave function (a;). By putting ex¬ 
pansion ([^ to the Hamiltonian Q one gets a general 
multi-band Hubbard-like Hamiltonian of the form 

^ = E E + E E cL^aia^ja^kasi. 

OL ij ijkl OLj3^5 

( 3 ) 

The parameters and are appropriate matrix 

elements of the original Hamiltonian Q 

= I dx (x) Ho will {x), (4a) 

= I / dxW["\x)Wf\x)wi^\x)wl^\x). 

(4b) 

The diagonal parameters represent average single¬ 
particle energies in the states described by Wannier func¬ 
tions wl°‘\x). In contrast, the off-diagonal terms t^^o 
are the tunneling amplitudes between lattice sites. Note, 
that due to the properties of Wannier functions the tun¬ 
neling between lattice sites does not change the orbital 
a. Since Wannier functions form a complete basis for the 
single-particle problem the Hamiltonian ^ is completely 
equivalent to the Hamiltonian 0. The advantage is that, 
in this picture, we are dealing with localized states and 
therefore it is a quite clear procedure leading to the sim¬ 
plified models. 

The standard BH model is obtained in the so-called 
tight-binding limit, i.e. when all inter-band and inter¬ 
site interactions can be neglected and only the nearest- 
neighbor tunneling is important. In the experiments 
with ultra-cold bosons this scenario is achieved for deep 
enough lattices and weak enough interactions when all 


FIG. 1: (Color online) Most relevant parameters of the ex¬ 
tended multi-orbital BH Hamiltonian <§ as functions of the 
depth of the optical lattice, (a) Ratio of the NNN tunneling 
amplitude t' to the tunneling amplitude t (opposite sign in¬ 
cluded). For a very shallow lattices the NNN tunneling is of 
the order of 10% of the tunneling t. (b) Ratio of the tunnel¬ 
ing amplitude t to the energy gap between bands A. Even 
for shallow lattices the tunneling amplitude is much smaller 
than the energy gap. (c) Ratio of the interaction-induced 
inter-site tunneling Uis to the local interaction parameter U. 
Due to the short-range character of the interactions all inter¬ 
site terms are negligible small, (d) Ratio of the inter-band 
coupling I/iB induced by interactions to the local interaction 
parameter U. In all plots, the parameters are calculated nu¬ 
merically from the exact shape of the Wannier functions. 


particles occupy the lowest Bloch band of periodic po¬ 
tential. In this limit only the on-site interaction in the 
lowest Bloch band C/™™ and the nearest-neighbor tun¬ 
neling are important. For simplicity we will use stan¬ 
dard notation U and t, respectively. Properties of the BH 
model are well known and were studied in many differ¬ 
ent contexts and with the help of different numerical and 
analytical methods HMg. 

Recently, due to the instantaneous experimental 
progress, many extensions of the standard BH model 
where introduced. Typically, extensions are needed 
whenever interactions between particles are enhanced 
and their long-range part starts to be important. In such 
scenarios the BH model is extended by an appropriate 
inter-site and inter-band terms. These path has opened 
very wide filed of research and has led to many interesting 
conclusions El- 

Here, we discuss a different regime of experimental pa¬ 
rameters where interactions between particles are still 
very weak but the periodic potential of the optical lattice 
is very shallow. Our aim is to derive the most relevant 
corrections to the standard BH model and find their con- 
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sequences in determining the position of the phase tran¬ 
sition point from the MI to the SF phase. 

The depth of the optical lattice is directly related to the 
energetical spreading of the Bloch bands (spatial spread¬ 
ing of the localized Wannier functions) and is responsi¬ 
ble for the narrowing of the energy gaps between them. 
Intuitively, for shallow lattices, one expects that the BH 
model needs to be extended by many different terms tak¬ 
ing these simple observations into account. In general, 
such a model is quite complicated and hard to analyze. 
However, subtle analysis shows that in the vicinity of the 
quantum phase transition point (for Filing p = 1) many 
of new terms can be completely neglected and the only 
correction which is relevant comes from the tunnelings to 
the NNN ^ 2 °^. In these particular case corrections that 
originate in the higher bands physics, as well as in the 
inter-site interactions, are less important and at the first 
approximation they can be omitted. 

To justify these non-obvious statements let us at first 
show that the influence of higher bands is negligble. 
The argumentation is based on the observation that 
the higher bands enter the game through interaction 
terms which lead to the promotion of particles from the 
ground to the higher bands. The most important term 
of this kind is related to the inter-band interaction Ui% = 
which is responsible for the inter-band transfer 
of two particles from the ground band |g) = (aQj)^|vac) 
to the first excited band |e) = (a|j)^|vac). In the sub¬ 
space of two particle states in given lattice site the local 
Hamiltonian (truncated to the two lowest bands) has the 
form 


Tiioc = 2A|e)(e| -P 2f7iB (|g)(e| + |e)(g|), (5) 


where A is the energy gap between lattice bands. When 
the inter-band interaction C/ib is taken into account the 
true local ground state |G) has some contribution from 
the excited state |e) and the squared projection 


l(e|G)P 


2a;2-p 2-t 


X = 


2C/ib 

A 


( 6 ) 


is a proper measure of the influence of higher bands. To 
estimate its value let us remind that, for one-dimensional 
case, the system undergoes the phase transition in the 
region where the ratio t/U > 0.25. This ratio is inde¬ 
pendent on the lattice depth and for given lattice can 
be controlled by adjustment of the interaction coupling 
constant g. The energy gap A between the ground and 
the excited band of the periodic potential highly depends 
on the lattice depth. The ratio t/A as a function of the 
lattice depth V is presented in Fig. [^. As one can note, 
the tunneling t is always much smaller than A and even 
for very shallow lattices it is not larger that 10% of A. At 
the same time the inter-band interaction (Jib is at least 
two times smaller than the on-site interaction [7 (Fig. 
W)- From these two facts and the following chain rule 

UlB _ UlB U t , , 

~ IT ' 1 ' A 


one hnds that in the vicinity of the phase transition 
point {U/t < 4) the interaction energy C/jb/A < 0.2 
for any lattice depth. Consequently, even for very shal¬ 
low lattices, the squared projection |(e|G)|2 < 4%. In 
addition, correction to the energy of the ground state, 
\6E\/A = 1 - y'l+x^ < 8%. 

This relatively small contribution from higher bands is 
caused mainly by an existence of a large energy gap A. 
In other words, all processes that lead to the excitation 
of particles to the higher bands are always off-resonant 
and they are effectively suppressed by the conservation 
of energy. Note however, that this argumentation breaks 
down for larger fillings due to the enhancement of the in¬ 
teractions, i.e. the interaction parameters are multiplied 
by the numbers of particles occupying appropriate states. 

The situation is quite different for processes acting 
within the ground-band of the lattice. In this case, the 
inter-site processes (tunnelings or interactions) couple 
states with relatively equal energies and therefore they 
are not suppressed by the conservation of energy. They 
lead to the non-local correlations that have a crucial im¬ 
portance in the vicinity of the transition point from MI 
to SF phase. The most relevant process of this kind is 
obviously the single-particle tunneling to the neighboring 
site controlled by t. To find the leading correction to this 
term for shallow lattices we should compare the influence 
of the NNN tunneling t' = with the influence of the 
interaction induced tunneling Uis = for f = j ± 1 

which transfers the particle to the neighboring site when 
the second particle is present nearby [8]. As it is seen 
in Fig. the inter-site amplitude Uis is always much 
smaller than U . It is worth noting that induced tunneling 
Uis is small directly due to the short-range character of 
mutual interactions. It is known that it can play a crucial 
role in the case of long-range interactions [5] . It can be 
shown that, in the vicinity of the phase transition, |t7is|/t 
is also essentially smaller than Therefore, for shal¬ 

low lattices, the most relevant correction comes from the 
NNN tunneling t' and we will study the model with this 
correction only. Obviously, further improvement of the 
model (which is beyond the scope of this article) would 
require other terms. One of the most important would 
be the tunneling induced interaction Uis- 

At this moment it should be emphasized that, due to 
the properties of Wannier functions, the tunneling am¬ 
plitude t' has an opposite sign when compared with tun¬ 
neling t. It means that the NNN tunneling effectively de¬ 
creases the kinetic energy of the particles. As explained 
later, this fact has significant and counterintuitive conse¬ 
quences on the stability of the insulating phase. 

In this way we finally obtain the extended BH model 
for shallow lattice of the form: 

U = {di-i -+- di+i) - 7 ' ^ al(0^-2 + di+2) 

i i 

+ ( 8 ) 
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Two tunneling amplitudes t and t' depend directly on 
the lattice depth V. In the tight-binding limit, when the 
NNN tunneling t' is neglected, the model is effectively 
controlled by the one parameter t/U. The value of this 
parameter, along with the density of particles p = IV/L, 
determines all properties of the ground-state of the sys¬ 
tem. From the experimental point of view this parameter 
can be controlled in two independent ways, i.e. by chang¬ 
ing the depth of the optical lattice or by changing the in¬ 
teraction coupling parameter g. As long as t' is excluded 
both experimental approaches are equally good from the 
model point of view. However, when the NNN tunnel¬ 
ing is taken into account the Hamiltonian is controlled 
by two independent parameters t/U and t'/U and in the 
second experimental scenario (when only interaction cou¬ 
pling g is changed) their ratio is fixed by the lattice depth, 
t' jt = const. Therefore, one can still study properties of 
the system as a function of the normalized tunneling t/U 
provided that the lattice depth V is fixed. In the first ex¬ 
perimental scenario when the lattices depth is changed, 
the ratio of both tunnelings t'/t varies (see Fig. This 
simple observation should be always taken into account 
whenever a dynamical quench through the phase transi¬ 
tion is considered. One should realize that two, seemingly 
equivalent, methods of controlling t/U can lead to differ¬ 
ent results due to the neglected NNN tunnelings. In the 
following studies, we assume that the parameters of the 
Hamiltonian are controlled for the given lattice depth via 
changing interaction strength g. 

At this point it is worth a reminder that the model 
is valid in the restricted range of parameters, i.e. for 
the first insulating phase (p = 1) and in the vicinity of 
the phase transition point. Nevertheless, the model is 
sufficient to determine the value of critical tunneling for 
different lattice depths. 


III. METHOD AND RESULTS 

As noted previously, here we assume that the lattice 
depth is fixed and the ratio t/U is tuned by changing the 
local interaction term U. In this way the phase diagram 
for different lattice depths is obtained and the critical 
value of the tunneling as a function of the lattice depth 
is determined. For a very deep lattice the ratio t'/t tends 
to 0 and one should expect to recover the limit of the 
standard BH model (tight-binding limit). 

The phase diagram of the model studied is obtained 
via a straightforward method based on the procedure 
of an exact diagonalization of the Hamiltonian EHIIII]. 
First, we fix the number of lattice sites L and the num¬ 
ber of particles N. We perform exact diagonalization of 
the Hamiltonian Q in full many-body Fock space (i.e., 
each lattice site can be occupied with 0,... ,N particles) 
applying periodic boundary conditions. In this way we 
find the ground state of the system and its energy 

E{N,L). Then, for commensurate filling p = N/L = 1, 
we calculate upper and lower chemical potential p± de- 


shallow lattice (V=3 Ep) tight-binding limit 


0.3 



• 

0.3 






• • 




• 

• 

^, 0.2 




^, 0.2 


• • 




• • 

• 




• • 


0.1 



• 

0.1 


• 

• 


0 0.1 0.2 0 0.1 0.2 


1/L 1/L 

FIG. 2: (Color online) The upper (p+) and lower (p-) lim¬ 
its of the insulating phase as a function of the inverse of the 
system size 1/L for shallow (left panel) and deep optical lat¬ 
tice (right panel) at particular point t/U = 0.25. The NNN 
tunneling t' is calculated directly from the shape of the Wan- 
nier functions. The solid lines are linear fits to the numerical 
data points. Extrapolation to the inhnite system size 1/L = 0 
gives the energy gap F of the insulating phase in the thermo¬ 
dynamic limit. Note that the energy gap F is larger for shal¬ 
low lattices. It means that, in contrary to the naive intuition, 
taking into account NNN tunnelings stabilizes the insulating 
phase. 

fined as: 

p+{L)=E{L + l,L)-E{L,L), (9a) 

p_{L)=E{L,L)-E{L-l,L). (9b) 

These two quantities strongly depend on the lattice size 
L and we are interested in their values in the thermo¬ 
dynamic limit of infinite lattice. Therefore, we diago¬ 
nalize the Hamiltonian for different sizes of the lattice 
L = 4,8,10,12 and we extrapolate data to the L oo. 
The extrapolation procedure is based on the observation 
that, for large enough lattice sizes, the chemical potential 
scales linearly with the inverse of the lattice size 1/L. In 
this way we obtain upper and lower bounds of the in¬ 
sulating lobe as a function of the normalized tunneling 
t/U. In this limit we also define the energy gap of the 
insulating phase F = p,+—p-. In Fig. [^we present an ex¬ 
ample data points and the extrapolated values of /i’s for 
t/U = 0.25 in the lattice of infinite depth (tight-binding 
limit, t'/t = 0) and with V = 3Efi, respectively. 

In principle, the quantum phase transition point from 
the MI to the SF phase occurs in the system when the 
energy gap F becomes equal to zero. Technically, this 
definition cannot be adopted directly to the numerical 
calculations due to the numerical uncertainty of F. Here, 
it is assumed that this uncertainty comes mainly from 
the extrapolating procedure to the thermodynamic limit. 
As a consequence, the position of the transition point tc 
is defined as the hopping amplitude at which the insu¬ 
lating gap becomes smaller than this uncertainty m- 
In this way the phase diagram for a given lattice depth 
can be obtained and the position of the quantum phase 
transition can be estimated. In Fig. the phase di¬ 
agrams for two extreme lattice depths V = SEn and 
U —^ oo (tight-binding limit) are presented. As one can 
note, the position of the tip of the insulating lobe is sub¬ 
stantially affected by the presence of the tunnelings to 
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FIG. 3: Shape of the first insulating lobe in the phase diagram 
of the system in the case of a very shallow lattice {V = SEr) 
and in the tight-binding limit when NNN tnnnelings can be 
completely neglected. Changing in the shape and the move¬ 
ment of the tip of the lobe is visible. The plot is obtained 
from the numerical data points by extrapolation to the ther¬ 
modynamic limit {L —>■ oo) as discussed in the text. 


Z) 


0.38 -o 
0.36 - o 
0.34 - 
’o.32 - 
0.3 - 
0.28 - 


tight-binding iimit 


10 


Lattice depth V (Er) 


12 


14 


0.11 

0.09 

^ 0.07 

0.05 

0.03 


FIG. 4: (Color online) Position of the qnantnm phase transi¬ 
tion point from the MI to the SF phase as a fnnction of the 
depth of the optical lattice. For shallow lattices the critical 
tnnneling tc is larger and the insulating phase is more stable 
(npper panel). In contrast to higher dimensions the criti¬ 
cal valne of the chemical potential decreases with the lattice 
depth (bottom panel). 
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the NNN. The position of the transition point estimated 
for different lattice depths is presented in Fig. For 
very strong optical potential the tight-binding limit is 
achieved. Moreover, it is interesting to note that, con¬ 
trary to a naive intuition, for more shallow lattices (when 
the NNN tunnelings are enhanced) the insulating lobe 


is enlarged. It comes from the fact that the tunneling 
t' has an opposite sign to the ordinary tunneling t and 
some kind of destructive interference of both processes is 
present in the system. This fact was already noticed for 
an analogous model of a two-dimensional system on the 
basis of the quantum rotor model |22j and in the mean- 
field approximation |23j . The position of the tip is also 
shifted in the direction of the chemical potential fj-c- As 
it is seen in Fig. for shallow lattices the critical point 
occurs for smaller chemical potential. This behavior is 
typical only for the one-dimensional case since for higher 
dimensions opposite result was predicted [5^ [53] . 


IV. CONCLUSIONS 


To conclude, in this paper the extended BH model de¬ 
scribing the system of ultra-cold bosons in a shallow one¬ 
dimensional optical lattice was studied. It was shown 
that including additional single-particle tunnelings to the 
NNN is sufficient to describe the properties of the system 
in the vicinity of the quantum phase transition point for 
low densities (first insulating lobe). For a given depth 
of the optical lattice, the NNN tunneling amplitude was 
calculated from the exact shape of Wannier functions. To 
find the critical behavior of the system, exact diagonal- 
ization of the Hamiltonian was used. For shallow lattices 
the first insulating lobe is enlarged due to the opposite 
sign of the NNN tunneling and the critical value of the 
chemical potential /ic is decreased. 

Results presented can shed some light on the problem 
of recent experiments and related theoretical works on 
the quench through the quantum phase transition point. 
Typically, such experiments are done by changing the op¬ 
tical lattice intensity and NNN tunnelings are neglected. 
In the first experiment with ®^Rb atoms confined in a 
three-dimensional optical lattice |3] the quantum phase 
transition was located around V ^ HEr. In such a case 
NNN tunnelings seems to be very small (see Fig. 
and in the first approximation can be neglected. Never¬ 
theless, it is worth noting that, whenever high precision 
measurement is done, these tunnelings introduce a quite 
big uncertainty of 2% on the value expected from the 
tight-binding limit (Fig. |^. This observation can have 
deep consequences for the problem of very precise quan¬ 
tum simulations performed in optical lattices [24) . One 
should remember that the results obtained can be sub¬ 
stantially affected by a presence of commonly neglected 
tunnelings to distant sites. As shown here, for high pre¬ 
cision measurements, the quench process should be done 
by tuning of the interaction constant g in a deep lattice 
with fixed intensity. 

The author thanks M. Gajda for fruitful comments 
and suggestions. This research was supported by the 
(Polish) National Science Center Grant No. DEC- 
2012/04/A/ST2/00090. 
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